Strength of Hydrogen Bond Network Takes Crucial Roles in the Dissociation Process of Inhibitors from the HIV-1 Protease Binding Pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions

Random walk with barriers: Diffusion restricted by permeable membranes

Restrictions to molecular motion by barriers (membranes) are ubiquitous in biological tissues, porous media and composite materials. A major challenge is to characterize the microstructure of a material or an organism nondestructively using a bulk transport measurement. Here we demonstrate how the long-range structural correlations introduced by permeable membranes give rise to distinct features of transport. We consider Brownian motion restricted by randomly placed and oriented permeable membranes and focus on the disorder-averaged diffusion propagator using a scattering approach. The renormalization group solution reveals a scaling behavior of the diffusion coefficient for large times, with a characteristically slow inverse square root time dependence. The predicted time dependence of the diffusion coefficient agrees well with Monte Carlo simulations in two dimensions. Our results can be used to identify permeable membranes as restrictions to transport in disordered materials and in biological tissues, and to quantify their permeability and surface area.Comment: 8 pages, 3 figures; origin of dispersion clarified, refs adde

Capturing the essence of folding and functions of biomolecules using Coarse-Grained Models

The distances over which biological molecules and their complexes can function range from a few nanometres, in the case of folded structures, to millimetres, for example during chromosome organization. Describing phenomena that cover such diverse length, and also time scales, requires models that capture the underlying physics for the particular length scale of interest. Theoretical ideas, in particular, concepts from polymer physics, have guided the development of coarse-grained models to study folding of DNA, RNA, and proteins. More recently, such models and their variants have been applied to the functions of biological nanomachines. Simulations using coarse-grained models are now poised to address a wide range of problems in biology.Comment: 37 pages, 8 figure

Mesoscopic model for DNA G-quadruplex unfolding

[EN] Genomes contain rare guanine-rich sequences capable of assembling into four-stranded helical structures, termed G-quadruplexes, with potential roles in gene regulation and chromosome stability. Their mechanical unfolding has only been reported to date by all-atom simulations, which cannot dissect the major physical interactions responsible for their cohesion. Here, we propose a mesoscopic model to describe both the mechanical and thermal stability of DNA G-quadruplexes, where each nucleotide of the structure, as well as each central cation located at the inner channel, is mapped onto a single bead. In this framework we are able to simulate loading rates similar to the experimental ones, which are not reachable in simulations with atomistic resolution. In this regard, we present single-molecule force-induced unfolding experiments by a high-resolution optical tweezers on a DNA telomeric sequence capable of adopting a G-quadruplex conformation. Fitting the parameters of the model to the experiments we find a correct prediction of the rupture-force kinetics and a good agreement with previous near equilibrium measurements. Since G-quadruplex unfolding dynamics is halfway in complexity between secondary nucleic acids and tertiary protein structures, our model entails a nanoscale paradigm for non-equilibrium processes in the cell.Work supported by the Spanish Ministry of Economy and Competitiveness (MINECO), grant No. FIS2014-55867, co-financed by FEDER funds. We also thank the support of the Aragon Government and Fondo Social Europeo to FENOL group. Work in J.R.A.-G. laboratory was supported by a grant from MINECO, No. MAT2015-71806-R).Bergues-Pupo, A.; Gutiérrez, I.; Arias-Gonzalez, JR.; Falo, F.; Fiasconaro, A. (2017). Mesoscopic model for DNA G-quadruplex unfolding. Scientific Reports. 7:1-13. https://doi.org/10.1038/s41598-017-10849-2S1137Arias-Gonzalez, J. R. Single-molecule portrait of DNA and RNA double helices. Integr. Biol. 6, 904 (2014).Burge, S., Parkinson, G. N., Hazel, P., Todd, A. K. & Neidle, S. Quadruplex DNA: sequence, topology and structure. Nucleic Acids Res. 34, 5402 (2006).Lam, E. Y., Beraldi, D., Tannahill, D. & Balasubramanian, S. G-quadruplex structures are stable and detectable in human genomic DNA. Nat. Commun. 4, 1796 (2013).Siddiqui-Jain, A., Grand, C. L., Bearss, D. J. & Hurley, L. H. Direct evidence for a G-quadruplex in a promoter region and its targeting with a small molecule to repress c-MYC transcription. Proc. Natl. Acad. Sci. USA 99, 11593 (2002).Endoh, T. & Sugimoto, N. Mechanical insights into ribosomal progression overcoming RNA G-quadruplex from periodical translation suppression in cells. Sci. Rep. 6, 1 (2016).Hänsel-Hertsch, R., Di Antonio, M. & Balasubramanian, S. DNA G-quadruplexes in the human genome: detection, functions and therapeutic potential. Nat. Rev. Mol. Cell Biol. 18, 279 (2017).de Messieres, M., Chang, J. C., Brawn-Cinani, B. & La Porta, A. Single-molecule study of G-quadruplex disruption using dynamic force spectroscopy. Phys. Rev. Lett. 109, 058101 (2012).Koirala, D. et al. A single-molecule platform for investigation of interactions between G-quadruplexes and small-molecule ligands. Nat. Chem. 3, 782 (2011).Long, X. et al. Mechanical unfolding of human telomere G-quadruplex DNA probed by integrated fluorescence and magnetic tweezers spectroscopy. Nucleic Acids Res. 41, 2746 (2013).Ghimire, C. et al. Direct Quantification of Loop Interaction and pi-pi Stacking for G-Quadruplex Stability at the Submolecular Level. J. Am. Chem. Soc. 136, 15544 (2014).Garavís, M. et al. Mechanical Unfolding of Long Human Telomeric RNA (TERRA). Chem. Commun. 49, 6397 (2013).Fonseca Guerra, C., Zijlstra, H., Paragi, G. & Bickelhaupt, F. M. Telomere structure and stability: covalency in hydrogen bonds, not resonance assistance, causes cooperativity in guanine quartets. Chemistry-A European Journal 17, 12612 (2011).Yurenko, Y. P., Novotn, J., Sklen, V. & Marek, R. Exploring non-covalent interactions in guanine-and xanthine-based model DNA quadruplex structures: a comprehensive quantum chemical approach. Phys. Chem. Chem. Phys. 16, 2072 (2014).Poudel, L. et al. Implication of the solvent effect, metal ions and topology in the electronic structure and hydrogen bonding of human telomeric G-quadruplex DNA. Phys. Chem. Chem. Phys. 18, 21573 (2016).Li, M. H., Luo, Q., Xue, X. G. & Li, Z. S. Toward a full structural characterization of G-quadruplex DNA in aqueous solution: Molecular dynamics simulations of four G-quadruplex molecules. J. Mol. Struct-Theochem. 952, 96 (2010).Islam, B. et al. Conformational dynamics of the human propeller telomeric DNA quadruplex on a microsecond time scale. Nucleic Acids Res. 41, 2723 (2013).Stadlbauer, P., Krepl, M., Cheatham, T. E., Koca, J. & Sponer, J. Structural dynamics of possible late-stage intermediates in folding of quadruplex DNA studied by molecular simulations. Nucleic Acids Res. 41, 7128 (2013).Li, H., Cao, E. & Gisler, T. Force-induced unfolding of human telomeric G-quadruplex: a steered molecular dynamics simulation study. Biochem. Bioph. Res. Co. 379, 70 (2009).Yang, C., Jang, S. & Pak, Y. Multiple stepwise pattern for potential of mean force in unfolding the thrombin binding aptamer in complex with Sr2+. J. Chem. Phys. 135, 225104 (2011).Bergues-Pupo, A. E., Arias-Gonzalez, J. R., Morón, M. C., Fiasconaro, A. & Falo, F. Role of the central cations in the mechanical unfolding of DNA and RNA G-quadruplexes. Nucleic Acids Res. 43, 7638 (2015).Linak, M. C., Tourdot, R. & Dorfman, K. D. Moving beyond Watson-Crick models of coarse grained DNA dynamics. J. Chem Phys. 135, 205102 (2011).Rebi, M., Mocci, F., Laaksonen, A. & Ulin, J. Multiscale simulations of human telomeric G-quadruplex DNA. J. Phys. Chem. B 119, 105 (2014).Stadlbauer, P. et al. Coarse-Grained Simulations Complemented by Atomistic Molecular Dynamics Provide New Insights into Folding and Unfolding of Human Telomeric G-Quadruplexes. J. Chem. Theory Comput. 12, 6077 (2016).Parkinson, G. N., Lee, M. P. & Neidle, S. Crystal structure of parallel quadruplexes from human telomeric DNA. Nature 417, 876 (2002).Bhattacharya, D., Arachchilageand, G. M. & Basu, S. Metal Cations in G-Quadruplex Folding and Stability. Frontiers in Chemistry 4, 38 (2016).de Lorenzo, S., Ribezzi-Crivellari, M., Arias-Gonzalez, J. R., Smith, S. B. & Ritort, F. A Temperature-Jump Optical Trap for Single-Molecule Manipulation. Biophys. J. 108, 2854 (2015).Smith, S. B., Cui, Y. & Bustamante, C. Optical-trap force transducer that operates by direct measurement of light momentum. Methods Enzymol. 361, 134 (2003).Mergny, J. L., Phan, A. T. & Lacroix, L. Following G-quartet formation by UV-spectroscopy. FEBS letters 435, 74 (1998).Torrie, G. M. & Valleau, J. P. Nonphysical sampling distributions in Monte Carlo free-energy estimation: umbrella sampling. J. Comput. Phys. 23, 187 (1977).Kumar, S., Bouzida, D., Swendsen, R. H., Kollman, P. A. & Rosenberg, J. M. The weighted histogram analysis method for free-energy calculations on biomolecules I. The method. J. Comput. Chem. 13, 1011 (1992).Evans, E. & Ritchie, K. Dynamic strength of molecular adhesion bonds. Biophys. J. 72, 1541 (1997).Dudko, O. K., Hummer, G. & Szabo, A. Intrinsic rates and activation free energies from single-molecule pulling experiments. Phys. Rev. Lett. 96, 108101 (2006).Friddle, R. W., Noy, A. & De Yoreo, J. J. Interpreting the widespread nonlinear force spectra of intermolecular bonds. Proc. Natl. Acad. Sci. 109, 13573 (2012)

Optical Trapping with High Forces Reveals Unexpected Behaviors of Prion Fibrils

Amyloid fibrils are important in diverse cellular functions, feature in many human diseases and have potential applications in nanotechnology. Here we describe methods that combine optical trapping and fluorescent imaging to characterize the forces that govern the integrity of amyloid fibrils formed by a yeast prion protein. A crucial advance was to use the self-templating properties of amyloidogenic proteins to tether prion fibrils, enabling their manipulation in the optical trap. At normal pulling forces the fibrils were impervious to disruption. At much higher forces (up to 250 pN), discontinuities occurred in force-extension traces before fibril rupture. Experiments with selective amyloid-disrupting agents and mutations demonstrated that such discontinuities were caused by the unfolding of individual subdomains. Thus, our results reveal unusually strong noncovalent intermolecular contacts that maintain fibril integrity even when individual monomers partially unfold and extend fibril length.National Institutes of Health (U.S.) (Grant GM025874)National Science Foundation (U.S.). CAREER (Award 0643745

Single-Molecule Force Spectroscopy: Experiments, Analysis, and Simulations

International audienceThe mechanical properties of cells and of subcellular components are important to obtain a mechanistic molecular understanding of biological processes. The quantification of mechanical resistance of cells and biomolecules using biophysical methods matured thanks to the development of nanotechnologies such as optical and magnetic tweezers, the biomembrane force probe and atomic force microscopy (AFM). The quantitative nature of force spectroscopy measurements has converted AFM into a valuable tool in biophysics. Force spectroscopy allows the determination of the forces required to unfold protein domains and to disrupt individual receptor/ligand bonds. Molecular simulation as a computational microscope allows investigation of similar biological processes with an atomistic detail. In this chapter, we first provide a step-by-step protocol of force spectroscopy including sample preparation, measurement and analysis of force spectroscopy using AFM and its interpretation in terms of available theories. Next, we present the background for molecular dynamics (MD) simulations focusing on steered molecular dynamics (SMD) and the importance of bridging of computational tools with experimental technique

Search for strongly interacting massive particles generating trackless jets in proton-proton collisions at s = 13 TeV

A search for dark matter in the form of strongly interacting massive particles (SIMPs) using the CMS detector at the LHC is presented. The SIMPs would be produced in pairs that manifest themselves as pairs of jets without tracks. The energy fraction of jets carried by charged particles is used as a key discriminator to suppress efficiently the large multijet background, and the remaining background is estimated directly from data. The search is performed using proton-proton collision data corresponding to an integrated luminosity of 16.1 fb - 1 , collected with the CMS detector in 2016. No significant excess of events is observed above the expected background. For the simplified dark matter model under consideration, SIMPs with masses up to 100 GeV are excluded and further sensitivity is explored towards higher masses

Search for dark matter produced in association with a Higgs boson decaying to a pair of bottom quarks in proton-proton collisions at root s=13TeV

A search for dark matter produced in association with a Higgs boson decaying to a pair of bottom quarks is performed in proton-proton collisions at a center-of-mass energy of 13 TeV collected with the CMS detector at the LHC. The analyzed data sample corresponds to an integrated luminosity of 35.9 fb(-1). The signal is characterized by a large missing transverse momentum recoiling against a bottom quark-antiquark system that has a large Lorentz boost. The number of events observed in the data is consistent with the standard model background prediction. Results are interpreted in terms of limits both on parameters of the type-2 two-Higgs doublet model extended by an additional light pseudoscalar boson a (2HDM+a) and on parameters of a baryonic Z simplified model. The 2HDM+a model is tested experimentally for the first time. For the baryonic Z model, the presented results constitute the most stringent constraints to date.Peer reviewe

Evidence for Top Quark Production in Nucleus-Nucleus Collisions

Peer reviewe

Observation of the B-s(0) -> X(3872)phi Decay

Using a data sample of proton-proton collisions at root s = 13 TeV, corresponding to an integrated luminosity of 140 fb(-1) collected by the CMS experiment in 2016-2018, the B-s(0) -> X(3872)phi decay is observed. Decays into J/psi pi(+)pi(-) and K+K- are used to reconstruct, respectively, the X(3872) and phi. The ratio of the product of branching fractions B[B-s(0) -> X(3872)phi]B[X(3872) -> J/psi pi(+)pi(-)] to the product B[B-s(0) ->psi(2S)phi]B[psi(2S) -> J/psi pi(+)pi(-)] is measured to be [2.21 +/- 0.29(stat) +/- 0.17(syst)]%. The ratio B[B-s(0) -> X(3872)phi]/B[B-0 -> X(3872)K-0] is found to be consistent with one, while the ratio B[B-s(0) -> X(3872)phi]/B[B+-> X(3872)K+] is two times smaller. This suggests a difference in the production dynamics of the X(3872) in B-0 and B(0)s meson decays compared to B+. The reported observation may shed new light on the nature of the X(3872) particle.Peer reviewe